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Abstract 

We consider the evolution of a system composed of TV non-interacting point particles of mass m in a 
container divided in two regions by a movable adiabatic wall (adiabatic piston). In this talk we discuss 
the thermodynamic limit where the area A of the container , the number N of particles, and the mass 
M of the piston go to infinity keeping and |y fixed. We show that in this limit the motion of the 
piston is deterministic. Introducing simplifying assumptions we discuss the approach to equilibrium and 
we illustrate the results with numerical simulations. The comparison with the case of a system with finite 
(A, N, M) will be presented. 
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1 Introduction 

The "adiabatic piston" is a well known problem in thermodynamics. Two gases enclosed in 
a cylinder are separated by an adiabatic movable piston. Initially the piston is at rest at the 
position Xq and the two gases on the left (— ) and on the right (+) of the piston are in thermal 
equilibrium with pressure and temperature T^. At time t = the piston is let free to move 
and the problem is to find the final equilibrium state, i.e. the position Xf of the piston and the 
state (p^Tf 1 ) of the two gases (fig.l.a). 
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figure 1: (a) The adiabatic piston, (b) The simple piston 

Experimentally the "adiabatic piston" has been used before 1940 [|J, to measure the ratio 7 = — of 
gases. From 1929 until today this ratio has been mesured p| by observing the damped "harmonic" 
motion of a simple piston under a constant external force (fig.l.b). The main result of the recent 
experiments is that one should distinguish between two regimes, corresponding to weak and strong 
damping with very different properties, e.g. 

• for weak damping the frequency of oscillations is close to the values computed assuming 
adiabatic oscillations. 

• for strong damping, and small enough amplitudes, the frequency correspond to isothermal 
oscillations. 

Since 1963 the adiabatic piston became a controversial conceptual problem because it led to the 
conclusion that the two laws of thermo statics are not sufficient to yield the final state [3] . In fact 
one is forced to solve the equations obtained from thermodynamics to find the final state [4]. 
In 1998 we started to investigate this problem using a very simple microscopical model [5 — 8]. 
The system consists of N point particles of mass m in a cylinder of length L and cross section 
A. The cylinder is divided into two compartments, containing respectively N- and N+ particles 
(N = + N + ), by an "adiabatic piston" of mass M S> m orthogonal to the wall 
with no internal degrees of freedom (fig.l.a). Initially the piston is fixed and the two gases are 
in thermal equilibrium. At time t = the piston is released and moves without friction, along the 
x-axis under the colisions with the gas particles. It is assumed that the gas particles move with 
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constant velocity, except for purely elastic collisions on the boundaries of the cylinder and on the 
piston. The system is isolated and thus the energy E is conserved. 

Since the components of the particles velocities perpendicular to the piston play no role in 
the dynamic of the piston we assume that all the particles have velocities parallel to the x-axis 
and all distribution functions are independent of the transverse coordinates. Therefore if v and 
V denote the velocities of a particle and the piston before a collision, then under the collision 

v-*i/ =2V -v + a{v-V) V ->V' = V + a(v-V) (1) 

where a =— (2) 

M + m v ; 

2 Numerical Results 

A systematic study of the time evolution of this model using numerical simulations have led 
to the conclusion that the evolution of the piston proceed essentially in two stages with totally 
different properties and time scales [5 — 8]. 

1. In a first stage caracterised by a time scale of order t\ = L \J~^ the evolution X(t) of the 
piston has the following properties (fig.2) 

(a) It is deterministic 

N± 

(b) Let R± = m ~j^j~- F° r fixed R±, the evolution is independant of M for M large enough 
(e.g. M > 10 4 m). Moreover 

• if R± < 4, the damping is weak. The frequency of oscillations and the damping 
coefficient increase with R±. 

• if R± > 10, the damping is strong. The evolution is independant of N± and M. 

(c) The system evolves adiabatically (i.e. without heat transfer through the piston) toward 
a state of "mechanical equilibrium", where the pressure are equal p + = p~ = Pf, but 
the temperature different. In particular the temperature of the gas increases under 
compression. 

2. In a second stage, caracterised by a time scale of order t% = — t\, which is very large if one 
considers realistic numbers, the evolution X(t) is caracterised by (fig.3) 

(a) It is stochastic , with pressures p^ fluctuating around pf. 

(b) It depends strongly on M, but not on N±. More precisely the following scaling relation 
is observed Xm{P) = X(t/M) where X(t) is independant of M and N±. 

(c) The system evolves with transfer of heat through the piston toward a state of " thermal 
equilibrium" where p + = p~ and T + = T~ . 

In particular in this stage the temperature of the gas decreases under compression. 

Let us remark that some time after "mechanical" or "thermal" equilibrium is reached, the piston 
starts oscillations which then disappear in time (fig.2). These oscillations reflect the fact that 
the velocity distributions of the gases are not maxwellian and they will be responsible for the 
approach to thermalisation with maxwellian distribution of velocities (see also [10] were similar 
results have been obtained). 

3 Thermodynamic limit 

Although the second stage is qualitatively understood and rather well described by Boltzmann 
Equation [8] (see fig3.b), the same is not true for the first stage. To gain insight into truly 
adiabatic evolution, one can either study the simple piston (fig.l.b) where no heat is involved [9], 
or consider the adiabatic piston in the thermodynamic limit [7] where A — > oo, N± — > oo, M — > oo 
with fixed values for 

L, R± = m-± and 7 = 2m—- (3) 

M M+m 
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In this section we shall presents this second approach to the study of adiabatic evolution. Let us 
first consider a finite system, i.e. A, N±, M are all finite. We start from the Liouville Equation on 
the complete phase space for N particles and the piston,with the dynamics defined by the elastic 
collisions (on the boundaries of the cylinder and the piston). Integrating Liouville equations over 
all variables except one we obtain the equations for the one-particle distributions of the left and 
right particles p (x,v;t) and of the piston ^(X,V;t). It involves the two-particle distribution 
function p± tP (x,v;X,V;t) for one particle and the piston: 

(d t + vd x )p~ (x,v;t)=6(x + L/2)vp~ (x,v;t)+I~ (4) 
(dt + vd x )p + (x,v;t) =-S(x- L/2)vp + {x,v;t) + /+ (5) 
d t 4>(V;t)=I (6) 



where 4>(V;t) = dXV(X,V;t) (7) 



I~(x,v;t) = dV (V -v)p- ;P (x,v';x,V';t) - dV (v - V)p- iP (x, v; x, V; t) 

J V J — oo 

/V f + OO 

dV (v - V)p +;P (x, v'; x, V'; t) - dV (V - v)p +jP (x, v; x, V; t) 

I(V; t) = A { f V dv(V- v) [p- surf {v\ V- 1) - pt urf (v, V; t)} 

r+oo 

+ dv (v — V) [pt urf W, V; t) - pt urf (v, V; t)] 



and 



-oo 



Psurf&V't) = / dxp± tP (x ) v;x,V;t). 

J — oo 



(8) 



In the thermodynamic limit a = is zero, and Eq.6 can be written as 

d t <P(V; t) = - 7 |- 1 ^ + °° dv (v - V) 2 pJ urf (v, V;t) — jT dv (v - V) 2 p+ urf {v, V; t) | . (9) 
With the initial conditions we have considered one can show that the solution is 



<t(X,V;t)=6(X-X(t))6(V-V(t)) (10) 

where {ivZLA [F - { v ;t) -F} ty.,t)] 

r+oo rV 

F 2 ~(V;t)= / dv(v-V) 2 pJ urf (v;t) F+(V;t)= / dv (v - Vfp^ rf (v;t) 

JV J —oo 

p ± , p (x, v; X, V; t) =p ± (x, v; t)V(X, V; t) (12) 
=p ± (X(t),v;t)S(V - V(t)) = pf urf (v;t)6(V - V(t)) (13) 

In conclusion in the thermodynamic limit M — > oo, the evolution of the piston is deterministic 
given by the coupled equations 4 — 5, 10 — 16. -Eg. 13 can be written in the form 

j t V = £ (p- - P + ) - (\ + (V) + \-(V)) V (14) 



where p represent the pressure exerted by the particles which are going to hit the piston and 
A ± (y) are strictly positive functions which represent the friction coefficients : 

dvpt urf (v,t)v 2 -2m J dvpj urf (v;t)v 2 (15) 

and 

X ± (V) = ±l m ^[F±(V) - F±(0)] (16) 
Similary introducing the energy per unit area of the left/right compartment, 

/+oo J f 

dv p^ (v; t) -mv 2 , /9 ± (t>;i) = / dxp ± (x,v;t) 
-oo ^ J 

we obtain 

^e ± = ±2mF±(V)V 

which shows that the change in the energy of the gases is entirely due to the work done by the 
piston: the evolution is therefore adiabatic in the thermodynamic limit. 



4 Infinite cylinder (L = oo) and thermodynamic limit (M = oo) 

If the cylinder is infinite and if the initial conditions are such that \p~ — p + \ is sufficiently 
small, we can expect that the recollisions of particles on the piston will be negligible. In that case 
the functions F 2 (V) will be independant of t. Therefore the piston will evolve monotonically to 
the stationary state with constant velocity V solution of 

r+oo _ rV 

_ dvp-{v){v -V) 2 = dvp + (v) (v -V) 2 (17) 

JV J-oo 

and the approach to the stationnary state is exponentially fast with time constant tq = j, where 
A = \ + (V = 0) + \~(V = 0). Moreover the stationnary velocity does depend on the value 
and the time necessary to reach this stationnary state is proportionnal to (tti^j) ■ 



5 Thermodynamic limit for the piston (M = oo), but cylinder 
with finite length (L < oo) 

In this section, to simplify the notation, we take N + = N_ = N and R = m^. In the case of 
a finite cylinder, recollisions of particles on the piston imply that P sur t(v; t) depend explicitly on 
time and we have thus to study the complete set of coupled equations. 

We remark first that if the friction coefficient A is sufficiently large the piston will reach almost 
immediately the constant velocity V given by (20) and this uniform motion should proceeed until 
the induced "shock wave" reaches the boundary. This expected result is well confirmed by the 
numerical simulations for R = larger than 4 (fig.2a) . 

Since we can not solve the complete set of coupled equations we restrict ourselves to a qualitative 
discussion. If p + = p~ then X(t) = Xq (recall that (M = oo)). Therefore if the initial conditions 
are such that \p + — p~\ is small, we expect that \X(t) — Xq\ ^ 1 and we introduce the following 
assumption 

_/V_ 1 N j_ 1 

pTurfU * (-f) W) = p-W; PL-fit) * (^) r3 ^y = p + ^) 

p~(t) = 2e-(t)-^ = p-(t)k B T-(t); p + (t) *L 2e+(t) - * = p + (t)k B T+(t) 
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Moreover since |V| is small we replace the friction coefficients X^iV) by A ± (0) computed from 
maxwellian distributions, i.e. 



X± = ml kBmp±V ^ (18) 
With these approximations we have two constants of the motion 

T~ + T + + j^V 2 = 2T (19) 

VWx - VT+(L - X) = C. (20) 
We can thus find the equilibrium point, which is independant of N and M: 

^j-V(L-x f y = ^c 

Tj = 2T ^- (21) 
T+ = 2T (1 - ^) 
Finally linearising the equations around the equilibrium point yield 

x = — WqX — Xx (22) 



x =^-^H + r--t^° (24) 



The evolution is thus weakly damped if 



2-7T 1 

with period r = - (25) 

It is interesting to remark that wq coincide with the frequency for adiabatic oscillations [2] 

6 Conclusions 

Numerical simulations have shown that strictly adiabatic evolution can be obtained only in 
the thermodynamic limit for the piston. The study of a simple microscopical model has shown 
that in this limit the evolution of the piston is deterministic, caracterised by two regimes, weak 
and strong damping depending on whether R = M ^ s is small or large. For weak damping the 
frequency of oscillations coincide with the values obtained in numerical simulations. However the 
predicted value for damping is larger than the simulations by several orders of magnitude. For 
strong damping the initial uniform motion coincide with the simulation with the correct value for 
the friction coefficient. 

Using simplifying assumptions we obtained equations for the final equilibrium state in agreement 
with the numerical simulations; moreover one can show that the entropy of both gases increase 
in time, in agreement with the second law for adiabatic evolution. 

Finally it should be stressed that the adiabatic evolution that we have discussed is strictly hamil- 
tonian, however in all cases we obtain a damping mechanism. This damping mechanism is still 
not well understood and the question of whether one can predict the final equilibrium state from 
fundamental principles remain an open problem. 
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Figure 2: Adiabatic Evolution 

a) Position X(t) for R = 0.1, 0.2, 4, M = 10 5 , 
R = 10, M = 30.000 

b) Pressure p ± {t) for R = 10, M = 30.000 

c) Temperature T ± (t) for = 10, M = 30.000 



Figure 3: Approach to thermal Equilibrium 

a) Position 

b) Scaling property (Position) 

c) Temperature for R = 10, M = 100 
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